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Abstract 

The theoretical description of heavy quarkonium production at high-energy pp 
and ep colliders is reviewed. Predictions based on non-relativistic QCD factori- 
sation are confronted with recent charmonium and bottomonium data from the 
Tevatron and HERA. Potential shortcomings of the present theoretical analy- 
ses are discussed, and the prospects for quarkonium physics at the upgraded 
Tevatron and HERA colliders and at the LHC are summarised. 
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1 Introduction 



The production of heavy quarkonium states at high-energy colhders has been the sub- 
ject of considerable interest during the past few years. A wealth of new experimental 
results has become available, some of which revealed dramatic shortcomings of earlier 
quarkonium models. In theory, progress has been made on the factorisation between 
the short-distance physics of heavy-quark creation and the long-distance physics of 
bound state formation. The colour- singlet model has been superseded by a con- 
sistent and rigorous framework, based on the use of non-relativistic QCD (NRQCD) 
0, an effective field theory that includes the so-called colour-octet mechanisms. On 
the other hand, the colour-evaporation model of the early days of quarkonium 
physics has been revived |7HlO . 



However, despite the recent theoretical and experimental developments the range 
of applicability of the different approaches is still subject to debate, as is the quan- 
titative verification of factorisation. Because the quarkonium mass is still not very 
large with respect to the QCD scale, in particular for the charmonium system, non- 
factorisable corrections [p!IHl3[| may not be suppressed enough, if the quarkonium is 
not part of an isolated jet, and the expansions in NRQCD may not converge very well. 
In this situation, a global analysis of various processes is mandatory to assess the im- 
portance of different quarkonium production mechanisms, as well as the limitations 
of a particular theoretical framework. 

Much of the recent interest and theoretical development in quarkonium physics 
has been stimulated by experimental data from high-energy pp and ep colliders. The 
purpose of this review is to summarise the theoretical progress in understanding heavy 
quarkonium production, and to confront the predictions with recent charmonium and 
bottomonium data from the Tevatron and HERA. Forthcoming experimental results 
and theoretical advances, for example in the calculation of higher-order corrections, 
will sharpen the present picture and allow more quantitative tests of the theory. 
Therefore, this review will mostly focus on the generic features of the various quarko- 
nium production processes and highhght the prospects for quarkonium physics at 
future collider experiments. Quarkonium production has also been studied in various 
other processes like e~^e~ annihilation, Z decays, hadronic collisions at fixed-target 
experiments and B decays, and several excellent reviews have been presented in the 
literature which cover these aspects, see e.g. ||T^-[18|. 



The NRQCD factorisation approach to quarkonium production is introduced in 
Section ^ and compared to earlier quarkonium models like the colour- singlet model 
and the colour-evaporation model. Section ^ summarises the NRQCD predictions 
for charmonium production and polarisation at hadron colliders. The predictions are 
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confronted with recent experimental data from the Tevatron, and potential short- 
comings of the theoretical analysis are discussed. Charmonium production at HERA 
is reviewed in Section |, and the prospects for quarkonium physics at the upgraded 
HERA collider are summarised. Section ^ concludes with an outlook on future work 
that will be required to test the theoretical framework more rigorously. 



2 Quarkonium Production Mechanism 

The production of heavy quarkonium at high-energy colliders involves two distinct 
scales. The creation of a heavy quark pair is a short- distance process on scales 
of the order l/mq or smaller and can be calculated in perturbation theory. The 
subsequent non-perturbative transition from the intermediate QQ pair to a physical 
quarkonium, on the other hand, involves long-distance scales of the order of the 
quarkonium size l/(mQt>) or larger. Provided that these two scales are well separated, 
l/{mQv) ^ l/rriQ, the formation of the bound state should be insensitive to the 
details of the heavy quark creation process, which is essentially local on the scale 
of the quarkonium size. It is thus intuitive to expect that long-distance and short- 
distance physics in quarkonium production can be separated such that binding effects 
factorise into universal non-perturbative parameters. 



2.1 Non-Relativistic QCD 

The factorisation approach to quarkonium production [§] provides a systematic frame- 
work for separating physics at long-distance and short-distance scales. It is based on 
the use of non-relativistic QCD (NRQCD) |jl9[, an effective field theory in which 
the heavy quark and antiquark are treated non-relativistically. A non-relativistic 
approach to quarkonium physics is motivated by the success of potential models in 
describing static properties of charmonia and bottomonia. The effective field theory 
framework of NRQCD allows a rigorous analysis of quarkonium physics which can be 
used to systematically improve earlier model approaches. 

The most general effective Lagrangian of NRQCD takes the form 

>Cnrqcd = i'^ (^Do + ^ + X^ (^Do - j X + ^Ught + 5C, (1) 

where ip and x ^'^'^ two-component spinor fields describing the heavy quark and an- 
tiquark, respectively, and £iight is the usual relativistic QCD Lagrangian for gluons 
and light quarks. The covariant derivative is = + ig^^-, where = {A^, A) 
is the SU(3) gauge field and g is the QCD coupling. The term 6C in (|1|) includes 
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all possible operators consistent with the symmetries of QCD and reproduces the 
relativistic effects of the full theory. The leading corrections in 5C are bilinear in the 
quark (or antiquark) field: 



5/:bmnear = ^ 1^^ 3^ ^ + ^ {D ■ qE - qE ■ D) 

H tb"^ (iD X gE — gE x iD) ■ a th -\ — ih"^ gB ■ aih 

Sm^Q 2mQ 

+ charge conjugate terms, (2) 

where E^ = and = ^e^^^G^^ are the chromoelectric and chromomagnetic 
components of the gluon field-strength tensor G^^ . The spin symmetry of the minimal 
NRQCD Lagrangian is broken by the two terms in (§) that contain the Pauli matrix 
a. The creation and annihilation of heavy quark pairs is taken into account by local 
four-fermion operators in 5C: 

5A-fermion = ~T {'^^ ^iX){X^ Ki') ■ (3) 

The factors /Cj, }C[ are products of a spin and a colour matrix and may also contain 
polynomials in the spatial derivate D. The operators in 5£4-femiion annihilate a QQ 
pair in a colour and angular momentum state specified by /C^ and create a QQ pair at 
the same space-time point in a state specified by /Cj. The dimensionless coefficients 
Cj, di in the effective Lagrangian are sensitive only to short distances of order l/mq or 
smaller and can be determined by matching scattering amplitudes in NRQCD with 
the corresponding amplitudes in full QCD. 

NRQCD calculations can, in principle, reproduce QCD results for quarkonium 
observables to any desired accuracy by adding correction terms to the effective La- 
grangian 6C and tuning the couplings to appropriate values. 



2.2 NRQCD factorisation 

Within the framework of NRQCD, the cross section for producing a quarkonium state 
H can be expressed as a sum of terms, each of which factors into a short- distance 
coefficient and a long-distance matrix element 0: 

da{H + X)=J2da{QQ[n]+X){0''[n]). (4) 

n 

The sum includes all colour and angular momentum states of the QQ pair, denoted 
collectively by n. The short-distance coefficients da are proportional to the cross sec- 
tions for producing a QQ pair in the state n and with small relative momentum. They 
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can be calculated perturbatively in the strong coupling as-f\ The non-perturbative 
transition probabilities from the QQ state n into the quarkonium H are given by 
vacuum expectation values of local four-fermion operators in NRQCD. Their generic 
form is 

{O''[n])^Y.(^\x^^nmW+X){H{X)+X\^p^]CU\0), (5) 
x,\ 

where the sum is over the quarkonium polarisation A and any number of light hadrons 
X in the final state. The factors /C„, K.'^ can contain products of colour and spin ma- 
trices as well as covariant derivatives. The colour and angular momentum quantum 
numbers n of the intermediate QQ pair need not equal those of the physical quarko- 
nium H. Soft gluons with energies of order rnqv or smaller can be emitted before 
the bound state is formed and change the colour and spin of the heavy quark pair. 
The effects of these soft gluons are included in the long-distance matrix elements 
(0^[n]).a 

The derivation of the factorisation formula (^) relies on the general space-time 
structure of heavy quark production to separate the short-distance effects of the hard 
momentum scale ~ tuq from the long-distance effects of the soft momentum scale 
~ Aqcd- This step, which is referred to as topological factorisation, uses standard 
factorisation methods of perturbative QCD |2^. The additional step required to 
derive the factorisation formula for the production of non-relativistic bound states 
involves separating the effect of the QQ relative momentum ~ mqv from the hard 
momentum scale rriQ. This can be achieved by expanding the hard scattering part of 
the amplitude as a Taylor series in terms of the heavy quark relative momentum. The 
long-distance factors generated by the Taylor expansion can be identified with the 
vacuum expectation values of local four-fermion operators in NRQCD. In the resulting 
expression all effects of the momentum scale rnqv or smaller are contained in 
the long-distance matrix elements {O^ [n]). They have to be calculated using non- 
perturbative methods or determined from experimental data. The short- distance 
coefficients da, on the other hand, involve only momentum scales of order ttlq or 
larger and therefore have a perturbative expansion in powers of as{mQ). A general 



prescription for calculating the short- distance coefficients da can be found in |21 



The factorisation of quarkonium cross sections has not been proven with com- 
plete rigour. Factorisation implies, for example, that soft gluons connecting the 
quarkonium and remnant jet parts of the scattering amplitude cancel for inclusive 
processes, but no explicit proof has been presented so far. Although it should be pos- 



^In the case of hadronic production, the short-distance cross sections have to be convoluted with 
parton distribution functions. 

^The dependence of the long-distance matrix elements and short-distance coefficients on the 
NRQCD factorisation scale has not been indicated and cancels in their product (^). 
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sible to extend the existing factorisation theorems of perturbative QCD to inclusive 
quarkonium processes, the apphcabihty of the NRQCD factorisation approach has 
not yet been estabhshed theoretically in all circumstances. 

The general factorisation formula contains an infinite number of terms. It 
provides a viable phenomenological framework only if the relative magnitude of the 
various non-perturbative transition probabilities can be determined systematically. 
NRQCD predicts a hierarchy among the matrix elements {O^ [n]) in terms of their 
dependence on the intrinsic heavy-quark velocity v. Using NRQCD power counting 
rules, the general expression (§) can be organised into an expansion in powers of v. 
To any given order in v only a finite number of terms contribute to the production 
cross section. 



2.3 Power counting 

The importance of the various long-distance matrix elements (O^ [n]) can be assessed 
using power counting rules. They depend upon the relative size of the three different 
low-energy scales in a non-relativistic bound state: the typical heavy quark three- 
momentum rnqv (in the quarkonium rest frame), the typical kinetic energy tuqv^ 
and Aqcd-S 



The standard NRQCD power counting (velocity scaling) rules have been 
derived from basic qualitative properties of non-relativistic bound states and from 
analysing the equations of motion for the NRQCD quantum field operators. For ex- 
ample, the number of heavy quarks in the quarkonium, / S'r tp'^ip ^ 1 implies ~ 
(rngv)^/^ since the quarkonium state can be localised within a region r ~ l/{mQv). 
Similarly, the kinetic energy operator / d^r ip"^ D"^/ (2mQ) ip has an expectation value of 
ttiqv'^, and so the spatial part of the covariant derivative has to scale like D ~ mgf . 
From the equation of motion for the vector potential one finds gE ~ ttLqV^ and 
gB ~ rnqv'^. According to these rules, the terms in 5£biiinear (0) are suppressed by 
0{v'^) relative to those of the leading NRQCD Lagrangian.0 

The power counting rules for the long-distance matrix elements {O^ [n]) in (^) 
can be derived by considering the Fock state decomposition of a quarkonium state 

''The power counting scheme for effective theories describing heavy-Hght mesons is more trans- 
parent. Because in a heavy-light meson the three-momentum and the energy of the heavy quark are 
both of order Aqcd , the expansion parameter is Aqcd /wq and the importance of the operators in 
the effective Lagrangian is determined by their dimensions. 

^The contributions from (S£4-fcrinion @) contain extra suppression factors Olaj.) from the short- 
distance coefficients di. 
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\H) in Coulomb gauge, 

\H) = V'Jq \QQ) + i^QQ, \QQg) + ■ ■ ■ (6) 

The dominant component \QQ) comprises a heavy quark pair in a colour-singlet state 
and with angular momentum quantum numbers "^^^^Lj that are consistent with the 
quantum numbers of the physical quarkonium. The higher Fock states, such as \QQg), 
contain dynamical gluons or light qq pairs. The heavy quark pair can be in either a 
colour-singlet or a colour-octet state with spin 5" = 0, 1 and angular momentum L — 
0, 1, 2, . . . All higher Fock states have probabihties suppressed by powers of v compared 
to \QQ). The \QQg) states with the highest probability of 0{v) are those that can be 
reached from the dominant \QQ) state through a chromoelectric interaction induced 
by the ip' {igA ■ d) ip term of the NRQCD Lagrangian. The chromoelectric transition 
satisfies the selection rules AL = ±1 and A^" = 0. Chromomagnetic transitions which 
proceed through the tp^ {gB ■ a) term in (5>Cbiiinear satisfy the selection rules AL = 
and AS — ±1. Higher Fock states \QQg) which can be reached from the dominant 
\QQ) state through a chromomagnetic interaction have probabilities of 0{v'^). Both 
chromoelectric and chromomagnetic transitions change the colour state of a QQ pair 
from colour-singlet to colour-octet, and from colour-octet to either colour-singlet or 
colour-octet. The chromoelectric transitions preserve the spin of the heavy quarks 
while chromomagnetic transitions are spin-symmetry breaking. 

The operator [n — 1{S), ^^'^^Lj] creates and annihilates a pointlike QQ pair 
in a colour-singlet (colour octet) state with angular momentum quantum numbers 
'^^'^^Lj. The magnitude of the corresponding vacuum expectation value {O^ [n]) can 
be determined from the power counting rules for the operator [n] and the proba- 
bilities of the quarkonium's Fock state components that have non-zero overlap with 
the quark pair in the state n. The general matrix element {O^ [n — 1(8), ^"^^^Lj]) 
scales as -y3+2L+2£;+4M^ where E and M are the minimum number of chromoelec- 
tric and chromomagnetic transitions required for the QQ pair to reach the state 
QQ [n = 1(8), ^'^"'"^Lj] from the dominant Fock state of the quarkonium H, and L is 
the orbital angular momentum quantum number. 

The standard NRQCD power counting rules have been derived assuming rnqv ^ 
mqv^ ~ Aqcd and as{mQv) ~ v. Such a hierarchy is likely to be realised for bottomo- 
nium states where niQ ~ 5 CeV and v"^ ~ 0.1. It is, however, not obvious a priori 
that it can be applied to the charmonium system where mq ~ 1.5 CeV and ^ 0.25, 
so that rnqv ~ 750 MeV. Depending on the relation between the low energy scales 
mQV,mQv'^ and Aqcdi alternative power counting schemes may be more appropriate 
for charmonium. If, for example, rnqv ~ Aqcd ^ rnqv'^ the typical momentum of dy- 
namical gluons would be ~ Aqcd, and chromomagnetic transitions may no longer be 
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suppressed with respect to chromoelectric transitions Alternative power count- 
ing rules can have important consequences for quarkonium phenomenology p3H25[|, 
as discussed in more detail in the subsequent sections. 



2.4 Colour-octet processes 

The NRQCD formalism implies that so-called colour-octet processes associated with 
higher Fock state components of the quarkonium wave function must contribute to 
the cross section. Heavy quark pairs that are produced at short distances in a colour- 
octet state can evolve into a physical quarkonium through radiation of soft gluons at 
late times in the production process, when the quark pair has already expanded to 
the quarkonium size. Such a possibility is ignored in the colour- singlet model 
where one assumes that only heavy quark pairs produced in the dominant Fock state 
form a physical quarkonium. 

The most profound theoretical evidence for the incompleteness of the colour- 
singlet model comes from the presence of infrared divergences in the production cross 
sections and decay rates of P-wave quarkonium. Consider for example the produc- 
tion of xj states. Within the colour- singlet model, xj is described by a colour-singlet 
QQ pair with angular momentum quantum numbers ^Pj. The cross section is as- 
sumed to factorise into a short- distance coefficient, which describes the production 
of a QQ [1, ^Pj] pair, and a single long-distance factor which contains all the non- 
perturbative dynamics of the bound state formation: 



The long-distance factor is related to the derivative of the radial wave function at the 
origin and can be determined from potential model calculations or from the xj decay 
width. At next-to-leading order, the short-distance cross sections for production (or 
decay) of P-wave quarkonia exhibit a logarithmic infrared divergence due to soft gluon 
emission, which can not be factored into |P'(0)p. In phenomenological applications 
of the colour-singlet model, an infrared cutoff has been introduced and identified with 
one of the low-energy scales present in non-relativistic bound states. It is evident, 
however, that the presence of infrared divergencies implies a failure of the simple 
factorisation assumption upon which the colour-singlet model is based. 

The NRQCD approach provides a natural solution to this problem [^. According 
to the NRQCD power counting rules, two terms contribute to the production cross 
section of P-wave states at leading order in the velocity expansion. For xj production 
one has 



dcr'^'^'^iXJ + X) = da{QQ [1, 'Pj] + X) (2J + 1) 



3Nc 



R'{0)\\ 



(7) 



daixj + X) 



da{QQ [l,3Pj]+X) {0^-'[i:'Pj]) 
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+ daiQQ[8,^S,]+X) {0>'-'[8,^S,]) + 0iv^). (8) 



The matrix element {O^-' [1, ^Pj]) can be related to the derivative of the radial wave 
function at the origin, and the first term in (j^) corresponds to the expression of the 
colour- singlet model. The second term represents a contribution to the cross section 
from a colour-octet mechanism. The short- distance factor is the cross section for 
producing a QQ pair in a colour-octet ^5*1 state, and the corresponding long-distance 
matrix element {O^-^ [8,^Si]) describes the probability for such a QQ [8,^5*1] pair to 
form a xj quarkonium. The colour-octet ^5*1 state can be reached from the dominant 
Xj Fock state through a single chromoelectric transition. According to the NRQCD 
power counting rules, the {O^' [1, ^Pj]) and {O^-^ [8, ^Si]) matrix elements scale with 
the same power of velocity ~ and must both be included in a consistent theoretical 
analysis. Other terms in the general factorisation formula are suppressed by relative 
order or more. In the NRQCD analysis, the infrared divergence in the colour- 
singlet short- distance cross section da{QQ [1, ^Pj] + X) at next-to-leading order is 
cancelled by a matching infrared singularity from the radiative corrections to the 
colour-octet matrix element {O'^-' [8,^5*1]). The inclusion of the colour-octet process 
proportional to {O^-' [8,^S'i]) is crucial to remove the infrared divergence from the 
production cross section.^ 

Colour-octet contributions are needed for a consistent description of P-wave quarko- 
nia, but they can be even more important phenomenologically for S-wave states like 
J/ip or ip{2S). According to the power counting rules, all colour-octet matrix ele- 
ments for the production (or decay) of S-wave quarkonia are suppressed by powers of 
the velocity compared to the colour- singlet contribution. The leading long-distance 
matrix element for J /if) production, for example, involves the operator O'^^'^ [l,^S'i], 
which annihilates and creates a cc pair with the same colour and angular momentum 
quantum numbers as in the dominant Fock state of the J / tp. From the power counting 
rules it follows that the corresponding vacuum expectation value (C'"^/'^ [1, ^^i]) scales 
as v^. According to the spin and angular momentum selection rules, the cc pair can 
reach the colour-octet states cc [8, ^Sq], cc [8, ^Si] and cc [8, ^Pj] through a single chro- 
momagnetic, a double chromoelectric, and a single chromoelectric transition, respec- 
tively. The corresponding vacuum expectation values {O'^^'^ [8,^So]), (C)-^/"^ [8,^5"!]) 
and (C-^/^ [8,^Pj]) all scale as v"^, i.e. they are suppressed by a power compared 
to the leading {O'^^'^ [1,^5*1]) matrix element. All other matrix elements {O'^l'^ [n]) 
scale as v^^ or smaller. In the non-relativistic limit f — > 0, the NRQCD description of 
S'-wave quarkonium production thus reduces to the colour- singlet model. However, 
as discussed in detail in the subsequent sections, colour-octet processes can become 

^For more comprehensive and pedagogical discussions of the treatment of infrared divergencies 
in production and decay of P-wave quarkonia see e.g. 
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significant, and even dominant, if the short- distance cross section for producing QQ 
in a colour-octet state is enhanced. 

NRQCD power counting rules predict a hierarchy among the various transition 
probabilities {O^ [n]). Truncating the expansion @ at a given order in the heavy- 
quark velocity, only a finite number of terms contribute to the production cross sec- 
tion. If the operator [n] creates and annihilates a QQ pair in a colour- singlet 
state with the angular momentum quantum numbers of the dominant quarkonium 
Fock state, its vacuum expectation value {O^ [n]) can be related to the quarkonium 
wave function at the origin and determined from the leptonic decay width or poten- 
tial model calculations. The number of independent non-perturbative matrix elements 
can be reduced further by using the approximate heavy-quark spin symmetry of the 
NRQCD Lagrangian 0. The calculation of the remaining NRQCD matrix elements 
that contribute at a given order in velocity expansion requires non-perturbative meth- 
ods. In the absence of lattice estimates^, the transition probabilities {O^ [n]) have to 
be treated as phenomenological parameters and determined from experimental data, 
see Section [31^ . 

Colour-octet processes are taken into account in the so-called colour- evaporation 
model where the cross section for a specific quarkonium state H is given as a 

universal fraction fij of the inclusive QQ production cross section integrated up to 
the open heavy quark threshold: 



The model assumes that any QQ pair, produced with small relative momentum, can 
form a physical quarkonium by emission of soft gluons. No constraints are imposed 
on the colour and angular momentum configuration. The cross section is hence dom- 
inated by QQ pairs produced at short distances with vanishing angular momentum 
and in either a colour-singlet or colour-octet state. 

The colour-evaporation model resembles some features of the NRQCD factorisa- 
tion approach, in particular due to the inclusion of colour-octet processes. However, 
the colour-evaporation model is not a rigorous theory, and there is no limit in which 
it could systematically reproduce a full QCD calculation. Moreover, it turns out 
that the assumption of a single universal long-distance factor is too restrictive from a 
phenomenological point of view. It implies, for example, a universal cr(xc)/o"(^/^) ra- 
tio, in confiict with experimental observations on charmonium production in hadron- 
hadron and photon-hadron collisions. 

"^First lattice results for quarkonium decay matrix elements have been presented in |^ . No lattice 
computations exist, however, for the non-perturbative matrix elements that describe quarkonium 





production. 
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3 Quarkonium Production at Hadron Colliders 



The first hadron colhder measurements of inclusive J/ipi '^{'^S) and Xc cross sec- 
tions from the UAl experiment |^| at the CERN SppS and from the CDF col- 



laboration at the Fermilab Tevatron could not separate charmonium produced 
directly in a hard-scattering reaction from charmonium produced in weak decays of 
5-hadrons. The early data showed some discrepancies when compared to the colour- 
singlet model, which predicts 5-decays to be the dominant source of charmonium 
at hadron colliders The substantial uncertainties in the theoretical calculation, 
however, made it difficult to draw definite conclusions |Q. A rigorous test of the 
colour- singlet model became possible when CDF presented a measurement of the di- 
rect ip{2S) cross section p4H36[|, where the contributions from S-decays had been 
removed using microvertex-detection. It turned out that only approximately 25% of 
i/^i^S) at the Tevatron come from the decay of i?-hadrons. Surprisingly, the direct 
iIj{2S) cross section at pti^l^) ^ 5 GeV proved to be orders of magnitudes larger than 



the leading-order colour-singlet model prediction P7|J5^ . A similarly striking dis- 



crepancy was observed in the direct J/ip cross section after contributions from both 
5-decays and radiative Xc decays had been removed.^ 

The cross section for quarkonium hadroproduction is given by 



da{pp ^ H + X) = Yl J dxidx2fi/p{xi)fj/p{x2] 



X 



^ da(t + j^QQ [n] + X) (O^ [n]), (10) 



where fi/p and fj/p denote the parton densities, and ij = gg, gq, gq, qq. In the colour- 
singlet model, the leading partonic subprocess g + g QQ [1, ^5*1] +g occurs at order 



37,38 



In retrospect, the shortcoming of the colour- singlet model to describe direct J/ip 
and ip{2S) hadroproduction can be understood by examining a typical leading-order 
Feynman diagram. Figure |T](a). At large transverse momentum, the two internal 
quark propagators are off-shell by ~ pf so that the parton differential cross section 
scales like da /dpi ~ l/pf , as indicated in the figure. On the other hand, when pt ^ 
2mc the quarkonium mass can be considered small and the inclusive charmonium cross 
section should scale like any other single-particle inclusive cross section ~ ^/Pt- The 
dominant production mechanism for charmonium at sufficiently large pt must thus be 
via fragmentation |^9[|, the production of a parton with large pt which subsequently 
decays into charmonium and other partons. A typical fragmentation contribution 

^Radiative x decays do not contribute to the ij{2S) cross section because the transition Xc 
ip{2S) is kinematically forbidden. 
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(a) leading-order colour-singlet: g + g ^ ccpS'j ] + g 




{b) colour-singlet fragmentation: g + g ^ [cc[^Si^^] + gg] + g 




(c) colour-octet fragmentation: g + g ^ ccpS'j ] + g 



(d) colour-octet t-channel gluon exchange: g + g ^ cc[^Sq^'' , "^Pj^^] + g 





Figure 1: Generic diagrams for J/ip and iIj{2S) production in hadron-hadron collisions 
through colour-singlet and colour-octet channels. 
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to colour- singlet J/ip production is shown in Figure |l](b). While the fragmentation 
contributions are of higher order in compared to the fusion process Figure [^(a), 
they are enhanced by a power p'l/{2m^^ at large pt and can thus overtake the fusion 
contribution at pt ^ 2mc. When colour-singlet fragmentation is included, the pt 
dependence of the theoretical prediction is in agreement with the Tevatron data but 
the normalisation is still underestimated by about an order of magnitude PO- 



indicating that an additional fragmentation contribution is still missing. It is now 
generally believed that gluon fragmentation into colour-octet ^5*1 charm quark pairs 
[^^ , as shown in Figure |l|(c), is the dominant source of Jjip and ip{2S) at large pt 
at the Tevatron. The probability of forming a, J/ip particle from a pointlike cc pair in 
a colour-octet ^5*1 state is given by the NRQCD matrix element (0^/^[8, ^^i]) which 
is suppressed by t>'^ relative to the non-perturbative factor of the leading colour- 
singlet term. However, this suppression is more than compensated by the gain in 
two powers of 0:^/77 in the short-distance cross section for producing colour-octet 
^Si charm quark pairs, as compared to colour-singlet fragmentation. At 0{v^) in 
the velocity expansion, two additional colour-octet channels have to be included, 
Figure |l|(d), which do not have a fragmentation interpretation at order but which 
become significant at moderate pt ~ 2mc The importance of the cc[8, ^5*0] and 

cc[8, ^Pj] contributions can not be estimated from naive power counting in as and v 
alone, but rather follows from the dominance of t-channel gluon exchange, forbidden 
in the leading-order colour-singlet cross section. The leading-order subchannels which 
contribute to J/ip and ipi2S) hadropro duct ion at order t>^ are thus 



9 + 9 
g + q/q 
q + q 



cc 
cc 
cc 



1,^5*1; 8,^5*1; 8,^5*0; 8, ^Pj + g, 
8,%; 8,%; 8,¥j] + g/g, 
8,^5*1; 8,^5*0; 8, ^Pj +g, 



(11) 
(12) 

(13) 



Explicit expressions for the parton cross sections can be found in [46| and [47 . 



3.1 Charmonium cross sections 

The different contributions to the J/ip transverse momentum distribution are com- 
pared to the CDF data |36| in Figure ^ As mentioned above, the colour- singlet model 



at lowest order in fails dramatically when confronted with the experimental results. 
When colour-singlet fragmentation is included, the prediction increases by more than 
an order of magnitude at large pt, but it still falls below the data by a factor of 
~ 30. The CDF results on J/ip production can be explained by including the lead- 
ing colour-octet contributions cc[8,^So], cc[8,^Si] and cc[8,^Pj], and adjusting the 
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corresponding non-perturbative parameters to fit the data.0 Numerically one finds 
the non-perturbative matrix elements to be of (9(10^^ GeV^), see Table |l|, consistent 
with the f ^ suppression expected from the velocity scaling rules. Similar conclusions 
can be drawn for 'ip{2S) production, see Figure |[ 



H 


m 






M3^5(8,%,3Po) 




1.16 GeV^ 


(1.19 ±0.14) X 10-2 


GeV^ 


(4.54 ± 1.11) X 10-2 GeV^ 


^(25) 


0.76 GeV^ 


(0.50 ±0.06) X 10-2 


GeV^ 


(1.89 ±0.46) X 10-2 GeV^ 


XcO 


0.11 GeV^ 


(0.31 ±0.04) X 10-2 


GeV^ 





Table 1: NRQCD matrix elements for charmonium production. The colour-singlet ma- 
trix elements are taken from the potential model calculation of |48,49|. The colour-octet 
matrix elements have been extracted from the CDF data |36], where M^(8, ^Sq, ^-Po) = 
(O^[8,^5'o]) + A; (0^[8, ^PoD/^T-c- "^^^ errors quoted are statistical only. Parameters: 
CTEQ5L parton distribution functions [^ ], renormalisation and factorisation scale /i = 

1.5 GeV. The Altarelli-Parisi evolution has been included for the 



p2 + 477^2 g^]^^ 

cc[8,^5i] fragmentation contribution. See |5l| for further details. 



Colour-octet processes need to be included in a theoretically consistent description 
of Xc production. According to Eq.(||) both cc[l,^Pj] and cc [8, ^5*1] intermediate 
states contribute to the Xc cross section at leading order in the velocity expansion. 
While the discrepancy between the colour-singlet prediction and the experimental 
results is less dramatic than in the case of S-wave charmonia, the inclusion 
of the ^Si colour-octet process (with {O'^" [S,^Si\) adjusted to the value hsted in 
Table 0) significantly improves the theoretical description. The comparison with 
recent data is presented in Figure ^. The measured ratio of production cross sections 
o"(Xc2)/o"(Xci) = 0.96 ± 0.27(stat) ± O.ll(sys) |5^ is also in good agreement with the 
NRQCD expectation of 1.1 ± 0.2 |5l. 

The analysis of the CDF cross section data alone, although very encouraging, does 
not provide a conclusive test of NRQCD factorisation because free parameters have 
to be introduced to fit the data. However, if factorisation holds the non-perturbative 

^The colour-octet ^S'q and ^Pj channels have a similar pt dependence, and the transverse momen- 
tum distribution is sensitive only to a linear combination of the [O'^^^ [8, ^^o]) and [O"^^^ [%,^Pj]) 
matrix elements. 
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Figure 2: Colour-singlet and colour-octet contributions to direct J/tp production in 
pp ^ J/iIj + X bX the Tevatron [^/s = 1.8 TeV, pseudorapidity cut |r/| < 0.6)) com- 



pared to experimental data from CDF \2>t\. Parameters: CTEQ5L parton distribution 



functions [^]; factorisation and renormalisation scale /i = yPt+ 4m2; rric = 1.5 GeV. The 
leading logarithms (a^ lnp^/(2mc)^)"' have been summed by solving the Altarelli-Parisi evo- 
lution equations for the gluon fragmentation function. NRQCD matrix elements as specified 
in Table |l|. 

matrix elements, Table |I|, are universal and can be used to make predictions for vari- 
ous processes and observables. The polarisation signature of J/ip and ip{2S) at large 
transverse momentum is one of the single most crucial tests of the NRQCD approach 
and will be discussed in some detail in Section |3.3| . Ultimately, a global analysis of dif- 
ferent production processes will be needed to assess the universality of the colour-octet 
contributions. Such a program, however, requires a careful discussion of the theoret- 
ical uncertainties associated with the determination of the non-perturbative matrix 
elements from experimental data. For charmonium hadropro duct ion the uncertainties 
are mainly associated with the calculation of the short-distance cross sections and, as 
argued in Section ^^2] , have not yet been reliably quantified. 

The J /ip and 'ip{2S) transverse momentum distributions at the Tevatron can also 
be described by the colour-evaporation model [|^-|^ or the related soft-colour inter- 
action model |jlO[- Both approaches allow colour-octet cc pairs produced in gluon 
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Figure 3: Same as Figure |2| for ip{2S) production. 
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Pt (GeV) 

Figure 4: Same as Figure |2| for Xc production. Experimental data from []52| . 



19 



fragmentation to hadronise into charmonium and thus resemble the NRQCD cross 
section prediction at large pt. The colour-evaporation model, however, assumes un- 
suppressed gluon emission from the cc pair during hadronisation, which randomises 
spin and colour, and consequently predicts unpolarised quarkonium. The NRQCD 
approach, in contrast, incorporates the spin-symmetry of QCD in the heavy quark 
limit, which implies transverse polarisation of direct J/ip and ip{2S) at large pt- The 



polarisation prediction of NRQCD will be discussed in detail in Section 3.3 



3.2 Determination of NRQCD matrix elements 

The numerical values for the non-perturbative matrix elements obtained from the fit 
to the Tevatron data are sensitive to any effects that modify the shape of the charmo- 
nmm pt distribution. This is true in particular for the (C»^[8, ^^o]) and (C^[8,3Po]) 
matrix elements which are determined from data at relatively small pt ^8 GeV. The 



theoretical uncertainties include the small- a; behaviour of the gluon distribution |51 



the evolution of the strong coupling |51|, as well as systematic effects inherent in 



NRQCD, such as the inaccurate treatment of the energy conservation in the hadroni- 



sation of the colour-octet cc pairs [£5| . Moreover, higher-order QCD corrections are 
expected to play an important role, as discussed in more detail below. 

The inclusion of higher-order QCD corrections is required to reduce the theoretical 
uncertainty and to allow a more precise prediction of the charmonium cross sections. 
Next-to-leading order (NLO) calculations for quarkonium production at hadron collid- 
ers are presently available only for total cross sections [^J57[] . Significant higher-order 
corrections to differential distributions are expected from the strong renormalisation 
and factorisation scale dependence of the leading-order results Moreover, the 

NLO colour-singlet cross section includes processes like g + g ^ QQ[^, '^Si] + g + g 
which are dominated by t-channel gluon exchange and scale as ~ a^(2mQ)^/p^. At 
Pt ^ 2mQ their contribution is enhanced with respect to the the leading-order cross 
section. Figure |l](a), which scales as ~ (2mQ)^/pf . This is born out by preliminary 
studies [^] which include part of the NLO hadroproduction cross section and by the 
complete calculation of NLO corrections to the related process of quarkonium pho- 
toproduction |5^. The NLO colour-singlet cross section may be comparable in size 



to the colour-octet ^5*0 and ^Pj processes, which scale as ~ a^f^(2mQ)^/pf (see Fig- 
ure |l](d)), and affect the determination of the corresponding NRQCD matrix elements 
from the Tevatron data. A full NLO analysis is however needed before quantitative 
conclusions can be drawn. 

Another source of potentially large higher-order corrections is the multiple emis- 
sion of soft or almost coUinear gluons from the initial state partons. These corrections. 
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as well as effects related to intrinsic transverse momentum, are expected to modify 
the shape of the transverse momentum distribution predominantly at relatively low 
values of pt <^ 2mQ. Initial state radiation can be partially summed in perturbation 
theory ||60|, but so far only total cross sections have been considered in the literature 
16T |. An estimate of the effect on the transverse momentum distribution should be 



provided by phenomenological models where a Gaussian fcf-smearing is added to the 
initial state partons. The result of these calculations, however, not only depends on 
the average (kf), which enters as a free parameter, but also on the details of how the 
smearing is implemented. Moreover, a lower cut-off has to be provided which regulates 
the divergences at pt = 0.Q Qualitatively, the inclusion of a Gaussian fcj-smearing 
leads to an enhancement of the short distance cross section at small pt, which results 
in smaller values for the fits of the {0^[8, ^Sq]) and {0^[8, ^Pq]) NRQCD matrix ele- 
ments [^2|,Q. The actual size of the effect, however, turns out to be different for the 
two models studied in the literature. 

An alternative approach to treat the effect of initial state radiation is by means 
of Monte Carlo event generators which include multiple gluon emission in the parton 
shower approximation. Comprehensive phenomenological analyses have been carried 
out for charmonium production at the Tevatron and at the LHC [|6^ -|65| using the 
event generator Pythia supplemented by the leading colour-octet processes ||63 



The inclusion of initial state radiation as implemented in Pythia leads to an enhance- 
ment of the short- distance cross section. The size of the effect is significantly larger 
than for the Gaussian fcf-smearing mentioned above, and it extents out to large pt- 
Consequently, the {0^[8,^So]) and {0^[8,^Po]) NRQCD matrix elements estimated 
from the Monte Carlo analysis of the Tevatron cross sections are significantly lower 
than those listed in Table ^ 

Recently, the impact of a non-vanishing initial state transverse momentum has 
been analysed in the context of the fcf-factorisation formalism. This approach has 
been developed to sum large logarithmic corrections in the high-energy limit s ^ 4m^. 
The summation generates initial state gluons with transverse momenta comparable to 
the hard scale mg. They are described by so-called unintegrated gluon distribution 
functions. The fc^-factorisation calculations ||6^j6^ indicate that the colour- singlet 
contribution to the direct J/ip cross section could be substantially larger than in the 
conventional coUinear factorisation approach. Moreover, the non-zero transverse mo- 
mentum of the initial state partons allows for the production of high-p^ colour-octet 
cc pairs without emission of hard gluons. The shape of the transverse momentum 

-'^"Using the NLO calculation for the total cross section ||5^, one can obtain the rough estimate that 
perturbative Sudakov effects should be confined below ^ 1 — 2 GeV for charmonium production 
at Tevatron energies. 
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distribution for the various colour-octet processes in the /cj-factorisation formahsm is 
thus different from the curves displayed in Figure The analyses suggest that 

the direct J/ip cross section might be dominated by the cc [8, ^Sq] and cc [8, ^Pj] chan- 
nels up to large transverse momenta of the order pt <^ 20 GeV. However, more work 
is needed to establish the quantitative predictions of the fct-factorisation approach in 
general. At present there is no reliable estimate of the corresponding theoretical un- 
certainties due to higher-order corrections and, in particular, due to the poorly known 
unintegrated gluon distribution. The impact of non-vanishing initial state transverse 
momentum has to be quantified in various other high-energy QCD processes before 
the importance of the different quarkonium production mechanisms can be assessed 
conclusively in this framework. 

To illustrate the uncertainty in extracting the NRQCD colour-octet matrix ele- 
ments from the J/ip, and Xc transverse momentum distributions at the Teva- 
tron, the results of various numerical analyses are collected in Tables |], |^ and ^. It is 
evident from the comparison that the theoretical error is not very well under control. 
The discrepancies between the analyses that rely on the leading-order coUinear fac- 
torisation calculation can be attributed to different choices for the input parameters 
and different prescriptions for the interpolation between the fusion cross section at 
low Pt and the fragmentation cross section including summation of ln(p^/(4m^) terms 
at large pt- More accurate data on the J/ip and ipi^S) cross sections could help to 
reduce some of the ambiguities. Ultimately, the calculation of next-to-leading order 
corrections to both colour-singlet and colour-octet processes is needed to improve the 
analysis and reduce the large factorisation and renormalisation scale uncertainty of 
the leading-order prediction. Effects beyond fixed order in perturbation theory, like 
the summation of soft gluon radiation, have to be addressed in a more rigorous and 
systematic way, and a partial summation of the NRQCD expansion might be needed 
to account for the inaccurate treatment of the kinematics in the hadronisation of the 
colour-octet cc pairs. 

Taking into account the large theoretical uncertainties of present analyses, most 
results appear to be consistent with the NRQCD power counting rules. The colour- 
octet matrix elements extracted from the fcj-factorisation calculation, however, do not 
follow the expected pattern. Such extreme values for the NRQCD matrix elements 
would indicate a serious flaw in our understanding of the power counting. 

It is instructive to compare the determination of the colour-octet matrix elements 
from the Tevatron data with analyses of other production processes. NRQCD ma- 
trix elements have been extracted from charmonium hadroproduction at fixed-target 
experiments and from i?-decays. The calculations include 0{as) corrections to the 
short-distance coefficients and, in the case of i?-decays, power corrections ~ AqcD/'^fe- 
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Reference 


PDF 


(0^/^[8,35i]) M//^(8,15o,=^Po) k 


LO collinear factorisation 


m 


MRS(DO) 


0.66 ±0.21 6.6 ±1.5 3 


HI 


CTEQ4L [P 
GRV-LO(94) 17| 
MRS(R2) PI 


1.06±0.14lJ;^^ 4.38 ± 1.15tJ:f^ 
1.12±0.14l[J|^ 3.90± 1.14^i:^^ 3.5 
1.40±0.22lJ:^^ 10.9 ± 2.07t?:^^ 


m 


MRST-LO(98) @ 
CTEQ5L [0 


0.44 ±0.07 8.7 ±0.9 

3.4 

0.39 ±0.07 6.6 ±0.7 


Parton shower radiation 


m 


CTEQ2L [0 
MRS(DO) 
GRV-HO(94) 


0.96 ±0.15 1.32 ±0.21 

0.68 ±0.16 1.32 ±0.21 3 

0.92 ±0.11 0.45 ±0.09 


[35 




CTEQ4M 


0.27 ±0.05 0.57 ±0.18 3.5 


/cj-smearing 


n 


CTEQ4M |71 


{h)[GeY] 
1 

1.5 


1.5 ±0.22 8.6 ±2.1 

3.5 

1.7±0.19 4.5 ±1.5 


m 


0.7 

MRS(D'_) [7C| 


1.35 ±0.30 8.46 ±1.41 

3 

1.5 ±0.29 7.05 ±1.17 


/c^-factorisation 


[37 




KMS 


[76| 


^0.04 ±0.01 ^6.5 ±0.5 5 



Table 2: NRQCD matrix elements for J/ip production obtained from various analyses 
of the J/V' transverse momentum distribution at the Tevatron. Values are given in units 
10"2 Qev3. (8,%,3Po) = {0^[8,^So]) + k {O^ [8,^ Pq]) /m^. The first error quoted 
is statistical, the second error, when listed, due to variation of the renormalisation and 
factorisation scales between = 1/2 Jpf + 4m^ and 2 Wpf + 4m^. 
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Reference 


PDF 




LO collinear factorisation 


46 


1 


MRS(DO) 17| 


0.46 ±0.1 1.77 ±0.57 3 


|51 


1 


CTEQ4L 
GRV-LO(94) 
MRS(R2) |7|] 


0.44 ± OMt^oil 1.80 ± O.bGf^oil 

0.46 ± 0.08lgi^ 1.60 ± 0.5ll[]i^ 3.5 

0.56 ± 0.1 ll[j:^^ 4.36 ± 0.96tJiJ 




1 


MRST-LO(98) 74| 
CTEQ5L g 


0.42 ±0.1 1.3 ±0.5 

3.4 

0.37 ±0.09 0.78 ±0.36 


Parton shower radiation 


B0] 


CTEQ2L |75l 
MRS(DO) l?! 
GRV-HO(94) 17| 


0.14 ±0.03 0.33 ±0.09 

0.11 ±0.03 0.28 ±0.07 3 

0.13 ±0.02 0.04 ±0.05 



Table 3: Same as Table ^ for ip{2S) production. Values are given in units 10-2 GeV^. 



Reference 


PDF 


(C»^=o[l,3Po]) [GeV^] (C»^=«[8,3Si]) [10-2 GeV=^] 


LO collinear factorisation 


46| 


MRS (DO) 


0.11 (input) 0.33 ±0.04 


|65| 


CTEQ4L |7T[ 


0.23 ±0.03 0.068 ±0.018 


|69| 


MRST-LO(98) ^ 
CTEQ5L |50| 


0.09 ± 0.01 (input) 0.23 ± 0.03 
0.09 ± 0.01 (input) 0.19 ±0.02 


/ct-factorisation 





KMS \7E\ 


0.11 (input) 0.03 ±0.01 



Table 4: Same as Table ^ for Xc production. 
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Representative values obtained for specific linear combinations of {0^[S,^So]) and 
(0'^[8,3Po]) are 



fixed-target hadroproduction < 



B decays 





^ 1.9 X 10-2 
^ 0.28 X 10- 


GeV^ 

2GeV^ U 




~ (i.5t;i) X 


10-2 GeV^ 




^ (0.5 ±0.5) 


X 10-2 GeV^ 


-'^-'3.4 


^ 2.4 X 10-2 


GeV^ 




^ 1.0 X 10-2 


GeV^ 
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Potentially large higher-twist effects at small pt may render the application of the 
NRQCD factorisation formula to quarkonium hadroproduction at fixed-target exper- 
iments unreliable. Nevertheless, the values obtained for the non-perturbative matrix 
elements are in qualitative agreement with the determination at the Tevatron. B- 
decays provide a more accurate testing ground for NRQCD factorisation. The lin- 
ear combination of ((9^[8, ^5*0]) and ((9^[8,^Po]) extracted from the branching ratios 
BR(i? tp + X) is consistent with Tevatron analyses that include soft-gluon effects 
in the parton-shower approximation. 

In subsequent numerical analysis the following ranges of values for the non-per- 
turbative matrix elements will be considered: 





(0^/'^[8,='5i]) 




= (0.3- 


2.0) X 10" 


-2 GeV^ ~ 






,iSo])+3.5(0^/'^[8 




= (1.0- 


10) X 10- 


-2 GeV^ 










= (0.1- 


1.0) X 10- 


-2 GeV=^ 


> (14) 




,iSo])+3.5(0^'(2^)[ 


8,3Po])/m2 


= (0.1- 


2.0) X 10" 


-2 GeV=^ 






(0>^-[8,35i]) 




= (0.1- 


0.5) X 10- 


-2 GeV^ 





Note that although the theoretical errors seem conservative, they do not accom- 
modate the non-perturbative matrix elements obtained from calculations in the kt 
factorisation approach. 



3.3 J/i/j and ij{2S) polarisation 

The analysis of J/ip and ip{2S) polarisation at large transverse momentum is one of 
the most decisive tests of the NRQCD factorisation approach. At large pt, the di- 

^^As alluded to above, the matrix elements obtained from calculations in the kt factorisation 
approach are not consistent with the NRQCD power counting rules. 
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rect J/ip and 'ip{2S) cross sections should be dominated by gluon fragmentation into 
colour-octet ^5*1 charm quark pairs, Figure [^(c). When pt ^ 2mc the fragmenting 
gluons are effectively on-shell and the intermediate heavy quark pairs transversely 
polarised. According to the standard NRQCD power counting rules, Section |2.3| , 
cc [8, ^5*1] evolves into J/tp or ip{2S) predominantly through a double chromoelectric 
transition which preserves the heavy quark spin. The dominance of gluon fragmenta- 
tion into colour-octet ^5*1 charm quark pairs at large pt and NRQCD spin-symmetry 
thus imply transversely polarised J/ip and ip{2S) at large transverse momentum jS^ . 

The polarisation can be measured through the angular distribution in the leptonic 
decay, 

oc 1 + a cos 9, (15) 



dcosO 

where 6 denotes the angle between the lepton three-momentum in the ip rest frame 
and the ip three-momentum in the laboratory frame. The polar-angle asymmetry a 
can take values in the range — 1 < a < 1, where a = 1 corresponds to transverse 
polarisation. 

Corrections to the asymptotic NRQCD prediction of pure transverse polarisa- 
tion can arise from spin-symmetry breaking chromomagnetic interactions in the non- 
perturbative evolution of the colour-octet ^5*1 charm quark pair. NRQCD power 
counting rules imply that spin-symmetry breaking transitions, which lead to random 
polarisation, are suppressed by a factor v"^. Hence their effect on a should not exceed 
~ 15% Additional corrections to transverse polarisation result from higher-order 
terms in the short-distance part of the gluon fragmentation function. After hard gluon 
radiation, the intermediate heavy quark pair can be in various spin and angular mo- 
mentum states and hadronise into unpolarised or longitudinal charmonium. The 



explicit calculation reveals that the effect is numerically quite small and that 
radiative 0{as) corrections can decrease the transverse polarisation by only about 
10%. The dominant source of depolarisation comes from the colour-octet ^Sq and 
^Pj fusion processes. Figure |l|(d). They are important for moderate transverse mo- 
mentum and need to be included to predict quarkonium polarisation in the pt-range 
probed by present data pT| , p^ . 

At 0{v'^) in the velocity expansion, the polarisation yield from all production chan- 
nels can be calculated unambiguously in NRQCD in terms of the non-perturbative 
matrix elements that have been determined from the unpolarised cross section.^ Fig- 
ure I displays the theoretical prediction for the asymmetry a in ip{2S) production 
at the Tevatron as a function of the ip{2S) transverse momentum. Both the 0{as) 
corrections to the fragmentation function and the (9(4m^/p^) corrections due to non- 

^^The calculation of the polarised cross section requires special care as interference contributions 
from intermediate cc pairs in different angular momentum states have to be included |8l],pl[ . 
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fragmentation contributions are included, but no attempt has been made to quantify 
spin-symmetry breaking effects. The shaded error band reflects the uncertainty in 
the determination of the NRQCD matrix elements from the unpolarised cross sec- 



tion. The matrix elements have been varied in the range specified in Equation 0, 
and the limiting cases that either [8,^So]) or [g^^Pj]) is set to zero 

in the combination M^^^\8, ^5*0, ^Pq) have been considered. No transverse polarisa- 
tion is expected at ~ 5 GeV, but the angular distribution is predicted to change 



drastically as pt increases. A first measurement from CDF does not support this 
prediction, but the experimental errors are too large to draw definite conclusions. 

The analysis of J/ ip polarisation is complicated by the fact that the experimental 
data sample ||83| includes J/ip that have not been produced directly but come from 
decays of Xc and ip{2S) mesons. The contribution from radiative x decays decreases 
the transverse J/tp polarisation at large pt-, but tends to be cancelled out by the 
contribution from ip{2S) decays []69| , |8^ . The overall polarisation yield for J /ip at large 
transverse momentum is thus not strongly diluted by decays of higher charmonium 
states. Again, the first experimental data do not confirm the NRQCD prediction [ |69| , 
see Figure ^. 

The absence of transversely polarized J /ip and ipi^S) at large transverse momen- 
tum, if confirmed with higher-statistics data at the Tevatron Run II, represents a 
serious challenge for the application of NRQCD to charmonium production. In the 
following, some uncertainties and potential shortcomings of the present theoretical 
analysis will be addressed: 

The uncertainty in the determination of NRQCD matrix elements {O"^ [8,^So]) , 
{O^ [8,^Si]) and {O"^ [8,^Pj]) translates into an uncertainty in the predicted yield 
of transversely polarized J/ip and ip{2S). In the large-p( limit, the ^Si production 
channel yields transverse polarisation, while the ^Sq and ^Pj channels both yield 
unpolarised quarkonia in this limit. A decreased value of the {O^ [8,^Si]) matrix 
element, relative to {O"^ [8, ^5*0, ^Pj]), would thus delay the onset of transverse polar- 
isation and reduce the discrepancy between NRQCD predictions and data. However, 
extreme values for the non-perturbative matrix elements outside the range specified 
in Equation |l^ would be in conflict with most of the analyses listed in Tables ^ and 
and moreover be inconsistent with conventional NRQCD power counting. 

Higher-order QCD effects in the short- distance cross section should be included to 
improve the accuracy of the theoretical prediction, in particular in the intermediate 
Pt region. Preliminary studies indicate large NLO effects for the unpolarized pt- 
distributions, see Section pTT| , but no attempt has been made so far to estimate the 
impact of higher-order corrections on the spin-dependence of the cross section. First 
studies within the fct-factorisation approach ||85| suggest that the J /ip and ip{2S) cross 
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Figure 5: Polar angle asymmetry a for ^{2S) production in pp ^{2S){-^ ^'^ jjL~) +X at 
the Tevatron as a function of pt compared to experimental data from CDF . Parameter 
specifications as in Figure ^. NLO corrections to the fragmentation contribution [^J5l]] have 
been included. The error band is obtained as a combination of the uncertainty (statistical 
only) in the extraction of the NRQCD matrix elements [Table |l| and the limiting cases that 
either (C'^[8, ^Sq]) or {0'^[S,^Pq\) is set to zero in the linear combination extracted from 
the data. 

section might not be dominated by gluon fragmentation but rather by the colour-octet 
^5*0 process up to pt ^ 20 GeV. Since the ^5*0 production channel yields unpolarised 
quarkonia, the overall J /ip and ip{2S) polarisation might be greatly reduced. However, 
as mentioned previously, the theoretical uncertainties in the fcj-factorisation prediction 
need yet to be quantified.]^ 

Heavy quark spin- symmetry is violated by higher-order terms in the NRQCD La- 
grangian, and longitudinal polarisation can arise if the binding of the charm quark pair 
into J /ip and ip{2S) proceeds through two chromomagnetic transitions. The suppres- 
sion of the spin-symmetry breaking chromomagnetic interactions follows from the con- 
ventional NRQCD power counting rules. These rules, however, have not been firmly 
established for the case of charmonium. Alternative power counting schemes [p!^ , p^ 
imply that the chromomagnetic transitions might not be suppressed with respect to 
chromoelectric transitions. As a consequence, the {O^ [8, ^Sq]) matrix element would 

^^The polarisation of Xcj has also been studied recently in the context of the /cj-factorisation 
approach Q. 
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Figure 6: Polar angle asymmetry a for prompt J /ip production (i.e. not from B-decay) in 
pp J/ip{-^ ii^n^) + X at the Tevatron as a function of pt compared to experimental 
data from CDF isl. From [lei. 



be enhanced relative to {O^ [8,^5*1]) and {O^ [8,^Pj]). Such a hierarchy is not in- 
consistent with the analyses of the unpolarised J/ip and ip{2S) cross sections, if one 
assumes that the linear combination M^(8, ^S'o, ^Po) is dominated by {O"^ [8,^So]). 
Since chromomagnetic transitions do not preserve heavy quark spin-symmetry the 
polarisation of J/ip and ip{'2,S) could be significantly reduced in the alternative power 
counting schemes. 

Finally, charmonium production mechanisms beyond NRQCD could be responsi- 
ble for J/ip and ipi'^S) depolarisation. For example, the leading- twist formalism of 
NRQCD factorisation does not include possible rescattering interactions between the 
intermediate heavy quark pair and a comoving colour field. It has been shown that 



such rescattering corrections could yield unpolarized quarkonium ||87| , |88[| . The anal- 
ysis |p7| , |88| , however, relies on several simplifying assumptions, and further work is 
needed to establish the importance of comover interactions for charmonium produc- 
tion at large Pf 

The distinctive NRQCD prediction of transverse J/ip and ip{2S) polarisation at 
large pt is not supported by first experimental data. However, the measurements 
still have rather large errors and no conclusions should be drawn until the statistics 
improves. If, in Run II of the Tevatron, the experimental results continue to dis- 
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agree with NRQCD predictions, one might have to conclude that either higher-order 
QCD effects could be crucial to describe the the spin-dependence of the charm cross 
section, or that alternative power counting schemes might be needed for the case of 
charmonium. 



3.4 Bottomonium production 

The application of NRQCD should be on safer grounds for the bottomonium system. 
As f ^ ~ 0.1 for bottomonium, higher-order terms in the velocity expansion (in partic- 
ular colour-octet contributions) are expected to be less significant than in the case of 
charmonium. Cross sections for the production of T and Xb states have been measured 
at the Tevatron in the region pt <, 20 GeV ||89H9ll • The leading-order colour-singlet 
model predictions underestimate the data, but the discrepancy is much less signifi- 
cant than in the case of charmonium. Given the large theoretical uncertainties in the 
cross section calculation, in particular at small pt ^ Mx, the need for colour-octet 
contributions is not yet as firmly established as for charmonium production. The 
inclusion of both next-to-leading order corrections and the summation of soft gluon 
radiation is required to obtain a realistic description of the T cross section in the 
Pf-range probed by present data. Such calculations have not yet been performed. 

The inclusive cross section for T(IS') production at the Tevatron is compared 
to the leading-order NRQCD calculation in Figure 0. The theoretical curves are 
based on the simple choice of non-perturbative parameters listed in Table ^ p^ . 
The colour-octet matrix elements have been determined from the CDF data for pt > 
8 GeV, where the contributions from direct T(IS') production and from radiative 
Xb decays are measured separately. Given the large uncertainty in the leading-order 
calculation of the short distance cross section, no attempt has been made to quantify 
the error of the non-perturbative parameters listed in Table |^. The ambiguities in the 
determination of the non-perturbative matrix elements for bottomonium production 
are best demonstrated by comparing the results in Table ^ with other numerical 
analyses |^^, collected in Table ^ for the case of direct T(IS') production. 



The impact of initial state gluon radiation on the T cross sections at the Tevatron 
has been estimated by adding a Gaussian /cf-smearing as discussed previously in 
the context of charmonium production. An average (kf) ~ 3 GeV and a ii"-factor 
~ 3 are found to bring the leading-order colour-singlet cross section in line with the 
inclusive T[1S, 2S) data [0. Similar results have been obtained within a Monte Carlo 
analysis |9^, leading to significantly lower fit values for the colour-octet NRQCD 
matrix elements than those determined from a leading-order calculation. Moreover, 
the Monte Carlo results imply that no feeddown from x states produced through 
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Figure 7: Colour-singlet and colour-octet contributions to inclusive T(15) production in 
pp T(15) -I- X at the Tevatron {^/s = 1.8 TeV, pseudorapidity cut |r/| < 0.4)) compared 



to experimental data from CDF |36|. Parameters: CTEQ5L parton distribution functions 



[pC|]; factorisation and renormalisation scale // = yPt +4m^; mb = 4.88 GeV. NRQCD 
matrix elements as specified in Table ^. 



colour-octet ^Si bb states is needed to describe the inclusive T cross section, in contrast 
to what is found at leading-order. The calculation of next-to-leading order corrections 
and a systematic treatment of soft gluon radiation within perturbation theory are 
required to resolve these issues. 

A preliminary measurement shows that inclusive T{1S) is produced unpolarised 
within the experimentally accessible range of transverse momentum. This result 
is consistent with NRQCD factorisation which predicts the onset of transverse T 
polarisation at Pt(T) ^ 15 GeV ||^, beyond what is probed by present data.Q 



^^Transverse polarisation of T states produced in p-Cu collisions at y/s = 38.8 GeV has been 
reported recently |^7| . The observation of a non-zero polarisation disproves the assumption that 
gluon emission from the bb pair during hadronisation randomises spin and colour, as implicit in 
the colour-evaporation model. On the other hand, the degree of transverse polarisation observed 
experimentally is smaller than predicted by a NRQCD calculation ||9^ , |99| . 
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Table 5: NRQCD matrix elements for bottomonium production. The colour-singlet ma- 
trix elements are taken from the potential model calculation of [^,49|. The colour-octet 

GeV iOtl, where 



matrix elements have been determined from the CDF data for pt > 
{0^\^,^ Sq\) = {0^\fi,^Po\) /m^ has been assumed for simplicity. Parameters: CTEQ5L 
parton distribution functions |5C], renormalisation and factorisation scale = ^Jpf + 4m^ 
and TUb = 4.88 GeV. 



Reference 


0^(i^)[8,35i]) (O^(i^)[8,i5o]) 




46 


1 


0.59 ±0.19 3.6 ±5 




93 


1 


2.0±4.l;[!J 13.6 ± 6.8t^°5^ 



Table 6: NRQCD colour-octet matrix elements for T(15') production obtained from various 
analyses of the transverse momentum distribution at the Tevatron. Values are given in 
units 10-2 GeV^. The first error quoted is statistical, the second error, when listed, due 
to variation of the renormalisation and factorisation scales between fi = 1/2 \Jpi + fn^ and 



3.5 Prospects for the Tevatron Run II and the LHC 

Run II at the Tevatron will provide a substantial increase in luminosity and will 
allow to determine the J/ip, V^(2S') and Xc cross sections more precisely and up to 
larger values of pt- A accurate measurement of the J/ip and ip{2S) polarisation at 
large transverse momentum will be the most crucial test of NRQCD factorisation. 
In addition, improved data on the J/ip and ^'(25') cross sections will help to reduce 
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some of the ambiguities in extracting the colour-octet matrix elements. 

With increased statistics it might be possible to access other charmonium states 
like the ric{nS) or the hc{nP). Heavy-quark spin symmetry provides approximate 
relations between the non-perturbative matrix elements that describe spin-singlet 
and spin-triplet states |0. The matrix elements for r]c{nS) are related to those for 
tp{nS), while the leading matrix elements for hc{nP) can be obtained from those 
for Xc{nP). Within NRQCD the rates for ri{nS) and h{nP) can thus be predicted 
unambiguously in terms of the non-perturbative matrix elements that describe the 
J/ip, ^(25) and Xc cross sections. Comparing the various charmonium production 
rates would hence provide a stringent test of NRQCD factorisation and heavy-quark 
spin symmetry. The cross sections for producing rjc and he at Run II of the Tevatron 
are large [[10CI| , |101|| , but the acceptances and efficiencies for observing the decay modes 
that can be triggered upon are in general small, and detailed experimental studies are 
needed to quantify the prospects. Other charmonium processes that have been studied 
in the literature include the production of D wave states |p.02|| , J/tp production in 
association with photons [ 103| , |104 | and double gluon fragmentation to J/ip pairs | 105| 



The larger statistics expected at Run II of the Tevatron will also allow to improve 
the measurement of the bottomonium cross sections. Yet undiscovered states like 
rih{lS) could be detected, for example in the decay rjh J/ip + J/ip |9^, and the 
associated production of T and electroweak bosons might be accessible ||106|| . If 



sufficient statistics can be accumulated, the onset of transverse T{nS) polarisation 
may be visible at Pt(T) ^ 15 GeV. 

A comprehensive test of the NRQCD factorisation approach in the bottomonium 
sector and a measurement of the T polarisation at pt ^ Mr will have to wait until 
the LHC starts operating. By that time, new experimental data from the Tevatron 
and HERA, as well as theoretical progress, e.g. in the calculation of higher-order 
corrections, will have significantly improved the present picture. Still, in order to 
provide some benchmark cross sections for quarkonium production at the LHC, the 
NRQCD predictions for selected charmonium and bottomonium states have been 
collected in Figure ^ Given the substantial uncertainty in the determination of 
the non-perturbative matrix elements from present data, the predictions should be 
regarded as order-of-magnitude estimates. 

The polarisation pattern predicted for direct T(IS') production at the LHC is 
presented in Figure ^, based on the NRQCD matrix elements of Table ^. Higher-order 
corrections to the gluon fragmentation function will lead to a small reduction of 
the transverse polarisation at large pt and should be included once data become 



-'^^For a collection of LHC cross section predictions obtained from a Monte Carlo event generator 
see |l07 . 



33 




10 20 30 40 50 

PxCGeV) 



Figure 8: Cross sections for charmonium and bottomonium production pp ^ H + X at 
the LHC {y/s = 14 TeV, pseudorapidity cut < 2.5). Parameters: CTEQ5L parton 
distribution functions |Q, factorisation and renormalisation scale n = \Jp'i + 4?tiq , mfic = 
1.5 GeV, mfe = 4.88 GeV. The leading logarithms (a^ Inp^ /(2mfe)^)"' have been summed 
by solving the Altarelli-Parisi evolution equations for the gluon fragmentation function. 
NRQCD matrix elements as in Tables |l],^. For bottomonium production the error band is 
obtained by varying the colour-octet matrix elements between half and twice their central 
value for illustration. 
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available. If the charmonium mass is indeed not large enough for a nonrelativistic 
expansion to be reliable, the onset of transverse T polarisation at pt ^ Mx niay 
become the most decisive test of the NRQCD factorisation approach. 
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Figure 9: Polar angle asymmetry a for direct T(IS') production in pp — > T{1S){-^ + 
X at the LHC as a function of pt- NRQCD matrix elements as specified in Table |5|, other 
parameters as in Figure 0. The error band reflects the limiting cases that either (0^[8, ^5*0]) 
or (C^[8, •^Pq]) is set to zero in the linear combination extracted from the data. 



4 Quarkonium production at ep colliders 

The analysis of charmonium and bottomonium cross sections at the high-energy ep 
collider HERA provides a powerful tool to assess the importance of the different 
quarkonium production mechanisms and to test the general picture developed in the 
context of NRQCD factorisation. Section [4.1| focusses on J/ip and 'ip{2S) photopro- 
duction, for which most of the experimental data have been collected. The impact of 
colour-octet processes and higher-order QCD corrections on cross sections and various 
differential distributions is reviewed, and the theoretical predictions are confronted 
with recent experimental results. Quarkonium production in deep-inelastic scattering 
is discussed in Section ^^2] . The high-statistics data expected at the upgraded HERA 
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collider will allow detailed studies of more exclusive observables like J/ip polarisa- 
tion, see Section [4.3| . Moreover, various other processes may be accessible, including 
Xc production, associated J/ip + production and the production of bottomonium 
states. The prospects for quarkonium physics at the upgraded HERA collider will be 



discussed in some detail in Section 4.4 



4.1 J/ip and ip{2S) photoproduction 

Quarkonium production in high-energy ep collisions at HERA is dominated by photo- 
production events where the electron is scattered by a small angle producing photons 
of almost zero virtuality. Experimental results have been reported for J/ip and ip{2S) 
photoproduction in a wide kinematical region, 30 GeV ^ ^J's^ ^ 200 GeV. The cross 
section is dominated by photon-gluon fusion, where the photon interacts as a point- 
like particle. Besides the direct production channel, charmonium photoproduction 
at HERA can also take place through resolved reactions, in which case the photon 
participates in the hard scattering through its parton content. Resolved processes 
contribute significantly to the lower endpoint of the charmonium energy spectrum 
and have been probed at HERA for the first time recently. Higher-twist phenom- 
ena like elastic/diffractive quarkonium production [p.08Hll0l are not included in the 
leading-twist calculation of NRQCD (see |^ for a detailed discussion) and have to 



be eliminated by either a cut in the J/ ip transverse momentum ^ 1 GeV, or by a 
cut in the J /ip energy variable z = pp ■ k^/pp ■ ^ 0.9 (in the proton rest frame, z 
is the ratio of the J/ip to j energy, z = E^/E^). 

Photoproduction mechanism The cross section for quarkonium photoproduction 
can be written as 



da{'yp^H + X)= ^ J dxidx2fi/jixi)fj/p{x2) 



«.ie{7.9.'?}' 

xY,da{i+j ^QQ[n]+X){0''[n]), (16) 

n 

where fi/-y and fj/p denote the parton distributions in the photon and the proton, 
respectively. For the direct photon processes i = 7 and /-^/^(x) = 5(1 — x). In the 
following, it is assumed that the J/ip transverse momentum is restricted to pt > 
1 GeV, in order to suppress the elastic contribution and higher-twist corrections in 
general. Generic Feynman diagrams which contribute to the parton cross section 



in the inelastic region away from pt = (or z = 1) are shown in Figure |T0|. At 
leading order in the velocity expansion, J/ip is produced through the colour- singlet 
channel |1|, Figure 0(a). Relativistic corrections at (9(f^) modify the large-^ and 



small-pt region |p.ll| - |114| | but can be neglected for pt I GeV. The next-to-leading 
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order QCD corrections to the direct photon process ||115| , p9| increase the cross section 
prediction significantly, in particular at pt ^ '^c as discussed in detail below. 

At Oiy'^) relative to the colour- singlet contribution, the J /ip can also be produced 
through intermediate colour-octet ^5*1, ^5*0 and ^Pj configurations ||116H118|| , see Fig- 
ure p!0|(b). Next-to- leading order corrections to the colour-octet channels are only 
known for the total photoproduction cross section (integrated over all z and pt) ||119| • 



(a) leading-order colour-singlet: 
direct 7: 7 + <? ^ ccpsf^'] -|- g 



resolved 7: + g ^ cc[^S[^^] + g 



7 WWV 




7 WV^^TOC- 




g v.Q.Q.O.O.Q.Q.Q.a 



(6) inelastic colour-octet: 



direct 7: 7-^5^ cc[^ S^^\^ P^^^] + g 



resolved 7: g^ + g ^ cc[^Si^^] + g 



WWVr 



vJ-0-O-P.O.' 






Figure 10: Generic diagrams for inelastic J/ip and ip{2S) photoproduction through colour- 
singlet and colour-octet channels. 



The importance of the colour-octet production processes follows from a kine- 
matical enhancement of the short- distance cross section, similar to the case of J/ip 
hadroproduction. The ^5*0 and ^Pj configurations can be produced through t-channel 
exchange of a gluon already at lowest order in a^. (For the octet channel this 
is true for the resolved process.) This leads to a potentially significant colour-octet 
contribution, in particular in the large- 2; region. 

The leading-order subchannels which contribute to inelastic J/ photoproduction 
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through direct photon processes are 

1 + 9 
-f + q/q 



cc 



cc 



1, ^5*1; 8,^5*1; 8/5*0; 8, ^Pj +9, 
8X; 8,%; 8;'Pj]+q/q, 



(17) 
fl8) 



where the initial-state parton originates from the target proton. For resolved photon 
processes, the subchannels are 



9 + 9 
g + q/q 
q + q 



cc 



cc 



cc 



Ij^'S'i; 8,^5*1; 8/50; 8,^Pj 
8,%; 8,%; 8,^Pj]+q/q 



+ 9, 



8,^5*1; 8/5*0; S/Pj 



+ 9, 



(19) 
(20) 
(21) 



where one of the initial-state partons originates from the photon and the other origi- 
nates from the proton. Explicit expressions for the parton differential cross sections 
are given in [^. The direct-photon processes dominate in the region z ^ 0.2, whereas 
resolved-photon processes become important in the region z <^ 0.2. At HERA en- 
ergies, photon-quark fusion can contribute about 10%-15% to the cross section at 
large z. Quark-gluon fusion constitutes about 20%-40% of the resolved cross sec- 
tion at 2 ^ 0.2 and becomes more important than gluon-gluon fusion at larger z. 
Quark- ant iquark fusion is always completely negligible. 

Fragmentation contributions in photon-gluon fusion | |120| - |122 ] exist only at the 
next order in and are suppressed for pt ^ 10 GeV. The cc[8, '^5'i] -channel is insignifi- 
cant for the direct-photon contribution, but dominates in resolved-photon interactions 
at ^ 5 GeV because it includes a gluon fragmentation component. The resolved 
photon amplitudes are identical to those relevant to J/ ip production in hadron-hadron 
collisions, and at HERA energies the relative importance of the various contributions 
as functions of pt is nearly the same as at Tevatron energies. 

The J/i/^ energy distribution The J/ip energy spectrum is a particularly sensitive 
probe of the different production mechanisms. The colour-singlet and colour-octet 
contributions to the differential cross section through direct and resolved photon 
processes are presented in Figure |lT] |p.l6|j47|| . The results are shown as a function of 
the scaling variable z = Ptp-Pp/p-^-Pp at a typical HERA energy of ^/s^p = 100 GeV in 
the inelastic region pt > 1 GeV. Adopting the NRQCD matrix elements as determined 
from the Tevatron J/ip data (see Section pl2| . Equation |1^, colour-octet processes are 
predicted to contribute significantly near the upper and lower endpoint of the energy 
spectrum. As mentioned previously, the colour-octet enhancement is a consequence 
of t-channel gluon exchange in the cc[8, ^5o] and cc[8, ^Pj] processes. The actual size 
of the effect, however, is difficult to quantify given the substantial uncertainties in 
the values of the (C»'^/'^[8, ^^0]) and {0-^/^[8,^Po\) matrix elements, as indicated by 
the shaded error band. 
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Figure 11: Colour-singlet (CS) and colour-octet (CO) contributions due to direct and re- 
solved photons to the J/ip energy distribution da/dz at the photon-proton centre-of-mass 
energy ^o^p 



/s^p = 100 GeV compared to experimental data from HI and ZEUS [124, 125 ], 
The shaded area represents the sum of all contributions with colour-octet matrix elements 
according to Equation |l^. The lines corresponding to separate colour-octet channels are 
plotted for {0-^/'('[8,%]) = (C'-^/'^[8, ^PoD/m,^ = 0.012 GeV^ The colour-singlet cross sec- 
tion is evaluated in leading order in a^. Other parameters: GRV(LO) proton and photon 
parton distribution functions |72,123|; = 200 MeV; factorisation and renormalisation 
scale fj, = 2m c- 



The theoretical predictions are compared to experimental results obtained by the 
HI pH and ZEUS |l25l collaborations at HERA in Figure |l|. The colour-octet 
enhancement in the region z ^ 0.6 is not supported by the experimental data. At 
first sight this could lead to rather stringent constraints on the colour-octet matrix 
elements (0"^/'^[8, ^S'o]) and (0'^/'^[8, ^Pq]) and might point towards an inconsistency 
with the values obtained for these matrix elements from other processes. However, 
the peaked shape of the 2;-distribution is derived neglecting the energy transfer of soft- 
parton emission in the non-perturbative transition of the colour-octet charm quark 
pairs. Such kinematic effects become important near the endpoint of the J / ip energy 
spectrum, and a summation of singular higher-order terms in the velocity expansion 
is necessary to predict the shape of the ^-distribution in the region z ^ 0.7 ||5^.[3 A 
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The situation is analogous to the well-known problem of extracting the CKM matrix element 
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quantitative study [|126|| reveals that the kinematic impact of soft-gluon emission in 
the hadronisation of colour-octet cc pairs at large z is not well under control, even 
after the summation of higher-order NRQCD corrections. A more stringent cut on 
the Jjil) transverse momentum is needed to remove the sensitivity of the prediction 
to processes which contribute in the elastic region at z ^ 1. At present, no reliable 
constraints on the size of colour-octet contributions in J photoproduction can be 
derived from the J jtl) energy spectrum near the upper endpoint. 

The low-z region is not expected to be sensitive to higher-order terms in the 
velocity expansion. Therefore, if the data could be extended to z ^ 0.3, an important 



resolved photon contribution should be visible [p.271J121| , |122|| , if the colour-octet matrix 
elements are not significantly smaller than extracted from the Tevatron data and 



expected from power counting rules. A preliminary measurement ||128|| indeed shows 
evidence for resolved processes, but the statistical error is still too large to establish 
the importance of colour-octet channels in the low-z region. 

NLO colour-singlet contributions The current experimental results on J/?/) 
photoproduction are adequately described by the colour-singlet channel alone, once 
higher-order QCD corrections are included [ 115 , ^ and once the normalisation uncer- 
tainty due to the choice for the charm quark mass and the QCD coupling is taken into 
account. The NLO calculation for the Jjil) energy distribution is confronted with the 
experimental data in Figure |12[ The shaded error band illustrates the large normali- 
sation uncertainty of the theoretical prediction: the upper boundary corresponds to 
a charm quark mass vnc = 1.3 GeV and as{Mz) = 0.1225, while the lower boundary 
has been obtained using rric = 1.5 GeV and as{Mz) = 0.1175.[] From Figure |l^ one 
can infer that the shape of the energy distribution in the colour- singlet channel is not 
very sensitive to the choice of input parameters. The theoretical prediction is in good 
agreement with the experimental results. A similar conclusion can be drawn from 
the analysis of the inelastic cross section as a function of the photon-proton energy. 
Figure 0. 

The next-to- leading order QCD corrections are crucial to describe the shape of 
the J/ip transverse momentum distribution. The NLO colour- singlet cross section 
includes processes like j + g ^ QQ[^, ^Si] + g + g which are dominated by t-channel 
gluon exchange and scale as ~ alrnq/p^. At pt ^ their contribution is enhanced 
with respect to the leading-order cross section, which scales as ~ a'^mq/p^. The 
comparison with the experimental data. Figure |13|, confirms the importance of the 
higher-order corrections. 



\Vub\ from the endpoint region of the lepton energy distribution in semileptonic B decay. 

^''The variation of the cross section with as has been studied in a consistent way by adopting the 
corresponding set of parton distribution functions [129|. 
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Figure 12: Next-to- leading order colour-singlet prediction for inelastic J/V' photoproduc- 
tion compared to experimental data from HI and ZEUS | 124|jl25 ,13C|. Upper Figure: J/^p 



energy distribution da/dz at the photon-proton centre-of-mass energy 



IP 



100 GeV. 



Lower Figure: inelastic J/ij: photoproduction cross section as a function of the photon- 



proton centre-of-mass energy J'sZ^. Parameters: MRST parton distribution functions |12 



factorisation and renormalisation scale /x = \f2m,c |59]. The error band is obtained by vary- 
ing the charm quark mass and in the range 1.3 GeV < < 1.5 GeV and 0.1225 
> as(Mz) > 0.1175, respectively. 
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Figure 13: Leading order and next-to- leading order colour-singlet prediction for the J /ip 
transverse momentum spectrum da /dpi at the photon-proton centre-of-mass energy ^/s' 



100 GeV compared to experimental data from HI and ZEUS [124,125]. Parameters as in 



Figure 12 



The colour-singlet channel provides a successful description of J/ip photoproduc- 
tion as far as present data on energy and transverse momentum distributions is con- 
cerned. The substantial normalisation uncertainties of the theoretical prediction, as 
well as the inaccurate kinematical treatment of soft gluon emission in the hadronisa- 
tion of colour-octet cc pairs, however, do not allow to place stringent constraints on 
the role of colour-octet processes. With the expected increase in luminosity at HERA, 
observables like J/ip polarisation and various other quarkonium final states will pro- 
vide more information on the importance of the different production mechanisms, see 
Sections 0| and 



t/?(2S') photoproduction While a non-relativistic treatment may be appropriate 
to describe J /if) photoproduction in many circumstances, significant relativistic cor- 
rections are expected for radial excitations like ^(25'). In the static limit f ^ 0, the 
a{ip{2S))/a{J/ip) cross section ratio is given by the ratio of the NRQCD long-distance 
factors (C^(25)^;l,35'^])/ {0-^/'('[l,^Si]) to all orders in perturbation theory. The 
colour-singlet matrix elements can be related to the quarkonium decay widths, so that 
aiipi2S))/a{J/ilj) = F^J25)/ j-j/^ ^ 0.41 ±0.05 Ip] in the static approximation. The 
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non-relativistic estimate is not confirmed by recent HERA measurements which yield 
a{ij{2S))/a{J/i)) = 0.210 ±0.058 IH and a{'ip{2S))/a{J/ij) = 0.242 ± 0.065 p^ , 



respectively. Higher-order terms in the velocity expansion will modify the theoret- 
ical prediction, their impact is, however, hard to quantify. Colour-octet processes 
might increase or decrease the cross section ratio, depending on the value of the non- 
perturbative matrix elements. No quantitative prediction is possible given the present 
theoretical error on {0^[8, ^Sq]) and {0^[8,^Pq]). Another source of relativistic cor- 
rections originates from NRQCD operators which involve total derivatives on fermion 



bilinears ||55| , [T6| . Although these operators are formally suppressed by powers of v"^ 
relative to the leading colour-singlet contribution, they are crucial to reproduce the 
hadronic phase space constraints of the quarkonium production cross section. Accord- 
ing to factorisation, all dependence on the specific quarkonium state resides within 
the NRQCD matrix elements. Since the short- distance coefficients are consequently 
expressed in terms of the heavy quark mass 2mQ rather than the quarkonium mass 
rriH, the physical phase space restrictions for producing a specific bound state are 
not manifest at fixed order in the velocity expansion. While 2mQ — rriH ~ Tnqv"^ 
is formally suppressed and numerically small for J/ip, the difference between parton 
and hadron kinematics can lead to a large ambiguity in predicting the cross section 
of radial excitations like ip{2S) where 2mc — ~ 1 GeV. The ambiguity can be 



resolved by summing certain higher-order terms in the velocity expansion The 
summation, however, leads to a non-perturbative distribution function, which can 
not be calculated except within models. A rough estimate [Q, obtained by assum- 



ing that the effective charm masses in the short-distance cross sections scale like the 
corresponding ip{2S) and J/ip masses, implies a significant phase space suppression 
factor of {mji^/m^Y ~ 0-6, consistent with the experimental observation. 

4.2 J/ifj production in deep-inelastic scattering 

Leptoproduction of quarkonium has been extensively studied in the context of the 
colour-singlet model ||133H136| and in the framework of the NRQCD factorisation 



approach including colour-octet processes [|137|J138|1 . At large photon virtualities Q^, 



the theoretical prediction is under better control than in the case of photoproduction 
because higher-order corrections and higher-twist contributions should become less 
important as increases. Moreover, diffractive processes [|110|] are expected to be 
suppressed at ^ 4mQ, and the total quarkonium leptoproduction cross section can 
be sensibly compared with the NRQCD prediction, contrary to photoproduction. The 
leading-order J/ip production process is 0{aas) and proceeds purely through inter- 
mediate ^5*0 and ^Pj colour-octet configurations. The leptoproduction cross section 
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in the high-energy hmit Q^,s ^ 4m^ ||137|| , 



lim, a{eP ^ e + J/ip + X) fdy fdx fg/p{x) 5{xys - Q^) 



mcQ"^ y \ m, 

{y is the momentum fraction of the J/ ip relative to the incoming electron) thus directly 
probes the NRQCD matrix elements (C-^/'^[8, %]) and (0-^/'^[8, ^Pq]) and provides a 
sensitive test of NRQCD factorisation. A more precise knowledge of {0'^^'^[8,^So\) 
and (0"^/'^[8, ^Pq]) would also help to clarify the role of colour-octet processes in J/ip 
photoproduction. At 0{aa1) J/ip is produced through "^Si colour-singlet and ^5*0, 
^Si and ^Pj colour-octet configurations. Results for the parton cross sections can 
be found in [|138|] . The 0{aa1) colour-octet processes are particularly important at 



large values of the inelasticity z. Note, however, that the upper endpoint region 
of the z-distribution in J/ip leptoproduction can not be predicted reliably without 
summing large higher-order corrections in the NRQCD velocity expansion, similar to 
phot opro duct ion . 

Experimental results on J/ip leptoproduction at HERA ||139|| are shown in Fig- 



ure |T^ compared to the leading-order NRQCD predictions. The theoretical cross 
sections include the 0{aa1) colour- singlet process and the 0{aas) colour-octet con- 
tributions. Taken at face value, the colour-singlet channel underestimates the data 
by a factor ~ 2 — 3. There is, however, a considerable amount of uncertainty in the 
cross section normalisation, in particular from the value of the charm quark mass, the 
strong coupling, and the choice of renormalisation and factorisation scales. On the 
other hand, adopting (0'^/^[8, %]) = 0.01 GeV^ and (0-^/^[8, ^Po])/^? = 0.005 GeV^ 
the sum of colour-singlet and colour-octet contributions exceeds the data by up to a 
factor of three. The shape of the distributions favour the prediction of the colour- 
singlet channel, but more experimental information extending to larger values is 
needed before firm conclusions can be drawn. 

4.3 J/ifj and iIj{2S) polarisation 

Polarisation measurements are powerful to assess the relative importance of colour- 
octet processes in different kinematic regions. First preliminary results [|140|] indicate 



a non-zero polarisation signature in J /if) photoproduction at HERA. With the ex- 
pected increase in luminosity, the analysis of decay angular distributions could pro- 
vide information on the relevance of the various quarkonium production processes, 
complementing the previously discussed polarisation measurement at hadron collid- 
ers. 
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Figure 14: Colour-singlet and colour-octet contributions to inelastic J ftp leptoproduction 
ep — > e J/V' X [137| compared to experimental data from HI | 139 |. The kinematic region is 
2 GeV^ < Q2 < 80 GeV^, 40 GeV < 1^ < 180 GeV. The variable y* denotes the rapidity 
of the J/^p in the 7*p centre-of-mass system. The invariant mass of the hadronic final state 
has been restricted to Mx > 10 GeV to suppress diffractive contributions. Parameters: 



GRV(LO) parton distribution functions |72,123]; factorisation and renormalisation scales 
fi = v^Q^ + 4m^; rric = 1.5 GeV. NRQCD matrix elements {0-^/'''[l,^Si]) = 1.1 GeV^ 



(0-^/^[8/5o]) = 0.01 GeV^, and {0-^/'^[8,^Po])/ml = 0.005 GeV^ 



The general decay angular distribution can be parametrised as 

— ^ ^ oc 1 + A cos^ 6^ + u sin 26 cos (f) -\ — sin^ 6 cos 2(j), 

ail 2 



(23) 



where 6 and refer to the polar and azimuthal angle of the /"^ three-momentum 
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with respect to a coordinate system defined in the J/ip rest frame. The parameters 
A, /X, V can be calculated within NRQCD as a function of the kinematic variables like 
z and pt Two examples are shown in Figure |T5|, the polar angle parameter 

A and the azimuthal angle parameter as a function of the energy fraction z and 
the transverse momentum pf, respectively. Each plot exhibits the result from the 
colour-singlet channel alone (dashed line) and the result after including colour-octet 
contributions. Since the decay angular distribution parameters involve ratios of cross 
sections, the dependence on parameters that affect the absolute normalisation, such as 
the charm quark mass, strong coupling, the renormalisation/factorisation scale and 
parton distribution, cancels to a large extent and does not constitute a significant 
uncertainty. 
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Figure 15: Colour-singlet and colour-octet contributions to the polar and azimuthal angu- 
lar parameters in the leptonic decay of J/ijj produced in photoproduction. (See [|^] for the 
definition of the Gottfried- Jackson and Collins-Soper reference frame.) The dashed line is 
the colour-singlet channel prediction. The shaded area bounded by the solid lines represents 
the sum of all contributions with colour-octet matrix elements according to Equation 



Parameters: GRV(LO) parton distribution functions |72,123]; factorisation and renormali- 
sation scale ^ = 2m^; rric = 1.5 GeV. 



The z-dependence of the polar angle asymmetry (Figure 0, left plot) is distinctive 
and can be used to discriminate between NRQCD and the colour-evaporation model, 
which always predicts unpolarised J/^.[^ The transverse momentum distribution, on 
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The predictions for the angular parameters in the endpoint region of the z distribution depend on 
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the other hand, can prove very useful to determine the relative magnitude of colour- 
singlet and colour-octet contributions: While colour-octet processes tend to predict 
unpolarised quarkonium in the pt region considered here, large and positive values 
of the azimuthal parameter u are expected for pt ^ 5 GeV if the cross section is 



dominated by the colour- singlet channel, see Figure |T5|, right plot. To observe the 
unique transverse polarisation signature of gluon fragmentation into colour-octet cc 
pairs one would have to isolate the resolved photon contribution at small 2; ^ 0.3 and 
measure the polar angular distribution for large transverse momenta pt ^ 10 GeV. A 
comprehensive analysis of the pt and z dependence of the angular parameters A, /x, u 
for different reference frames can be found in |^ . 



Polarisation in J/ip photoproduction through the colour- singlet channel has also 
been studied in the context of the fcj-factorisation approach ||141|| . The virtuality of 
the initial state gluon increases the fraction of longitudinally polarised J/ip, in par- 
ticular at large transverse momentum. The effect becomes significant at ^ 3 GeV, 
and could be probed at the upgraded HERA collider by analysing the polar angle 
asymmetry. The impact of colour-octet processes on the polarisation prediction in 
the /cf-factorisation approach has not been estimated so far. 

Information on the quarkonium production mechanism can also be obtained from 
the angular distribution in inelastic J/ip leptoproduction [ |135|J137|J138|| . The theoret- 
ical prediction at leading order only depends on the ratio of the NRQCD colour-octet 
matrix elements {0'^^^[8,^So\) and ((9"^/'^[8,^Po])- Longitudinally polarised J/ip is 
expected if the ^Pj channel dominates, while the ^5*0 channel leads to unpolarised 
quarkonium. The experimental result, A = 0.77lo;38 ||139|| , is consistent with colour- 
singlet dominance, which implies A ~ 0.5 [|139|| . The experimental uncertainties are, 
however, still too large to put stringent limits on the size of colour-octet contributions. 



4.4 Prospects for the upgraded HERA collider 

With the HERA luminosity upgrade, a wealth of new quarkonium data will become 
available. The existing measurements can be extended into new kinematic regions, 
and observables like J/ip polarisation and various other quarkonium final states will 
be accessible. The future analyses of quarkonium production at HERA offer unique 
possibilities to test the theoretical framework of NRQCD factorisation and to clar- 
ify the role of colour-octet processes. Some of the most interesting and promising 
reactions will be discussed in the following. 

the non-perturbative distribution functions that resuh from a summation of higher-oder corrections 
in the velocity expansion. Note, however, that the dependence is significant only if the angular 
parameter varies strongly in the endpoint region and disappears entirely if its distribution is flat . 
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Resolved photon processes As emphasised previously, the low-2; region of the 
J / ip energy spectrum in photoproduction is dominated by resolved photon processes, 
which resemble the J/ip production mechanisms relevant at hadron colliders. Away 
from the upper endpoint region, the energy spectrum is not expected to be sensitive 
to higher-order terms in the velocity expansion, and a measurement of the resolved 
cross section can provide a sensible test of the importance of colour-octet processes. 

Colour-singlet and colour-octet contributions to the low-z region of the energy 
spectrum are displayed in Figure left, at a photon-proton energy of W^p = 
170 GeV. Note that, because only a small fraction of the incoming photon energy 
participates in the hard interaction, the typical W^p values for resolved photon events 
within the detector acceptance are higher than for the medium- 2; analyses ||128|| . The 
right plot of Figure |16| shows the uncertainty of the leading-order colour-singlet cross 
section due to the variation of the charm quark mass and the renormalisation scale. 
Next-to-leading order corrections to inelastic J/ip photoproduction through resolved 
photons are not yet available. 

A preliminary measurement |p.28|| clearly establishes the existence of resolved pro- 
cesses at small z, but experimental errors and theoretical normalisation uncertainties 
do not allow to constrain the size of colour-octet contributions. To confirm the sig- 
nificant excess over the colour- singlet process observed at the Tevatron it will be 
necessary to analyse resolved photon processes at larger values of pt- 

Xc photoproduction The measurement of inelastic Xc photoproduction is a partic- 
ularly powerful way to discriminate between NRQCD and the colour-evaporation 
model. The assumption of a single universal long-distance factor in the colour- 
evaporation model implies a universal (T{xc)/cr{J/4') ratio, and a large Xc cross 
section is predicted in photon-proton collisions, similar to hadron colliders where 
(t(xc)/(x(J/V^) ^ 0.5 ig. In NRQCD, on the other hand, the cr(xc)/cr(J/^) ratio 
is process-dependent and strongly suppressed in photoproduction. While Xc can be 
produced copiously in hadron-hadron collisions through gg, gq and qq interactions, 
the direct photon process 7 + — > cc[l, ^Pj] + g is forbidden at leading order in ag 
due to the specific colour structure and charge conjugation invariance. Apart from 
higher-order corrections and resolved photon interactions, Xc can thus only be pro- 
duced through the colour-octet process 7 + (7 — > cc[8, ^^i] + g ||142| , |127|| . According to 
factorisation, and up to corrections of the cross section ratio a{xc)/o'{J/4') 

in direct photon interactions can be expressed in terms of the corresponding long- 
distance factors 

a(7P ^ XjX) a{^p ^ cc[8,^S^]X) {0^■'[8,^S^]) 



a{jp^J/^X) a{-fp ^ cc[l,^Si]X){OJ/'^[l,^Si]) 
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Figure 16: Direct and resolved photon contributions to the J/ip energy distribution da/dz 
at the photon-proton centre-of-mass energy ^s^yp = 170 GeV. Left: Colour-singlet (CS) and 
colour-octet (CO) contributions. The shaded area bounded by the solid lines represents 
the sum of all contributions with colour-octet matrix elements according to Equation 
The separate colour-octet channel is plotted for {0^/'f'[8,^So]) = {0-^/'l'[8,^Po]) /ml = 
0.012 GeV^. Other parameters as in Figure |l^. Right: Colour-singlet contribution only. 
The error band is obtained by varying the charm quark mass and the renormalisation scale 
in the range 1.3 GeV < nic < 1.5 GeV and rric < /^_r, < 2mc, respectively. 



(2J + 1)0.005, 



(24) 



with NRQCD matrix elements as listed in Table ||. A search for Xc production at 
HERA resulting in a cross section measurement or upper cross section limit would 
directly probe the colour-octet matrix element {0^-'[8,^Si]) and test the assumption 
of a single universal long-distance factor implicit in the colour-evaporation model. 

Photoproduction of spin-singlet states The inclusion of colour-octet processes 
is crucial to describe T]c{nS) and hc{nP) photoproduction. With regard to the P-wave 
state he, the colour-octet contribution is required to cancel the infrared divergence 
present in the colour-singlet cross section [ p.43|| . The production of rjc, on the other 
hand, is dominated by colour-octet processes since the colour-singlet cross section 
vanishes at leading-order due to charge conjugation invariance ||144| , |145|| , similar to 
Xc photoproduction. As discussed previously in the context of future quarkonium 
physics at the Tevatron, Section |3.5|, heavy-quark spin-symmetry allows to predict 
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the production of spin-singlet states like r]c{nS) and hc{nP) in terms of the NRQCD 
matrix elements that describe J/tp and Xc production. The cross sections for ric{nS) 
and hc{nP) photoproduction are sizeable [ p.43|jl44| |, but it is not obvious that these 
particles can be detected experimentally, even with high-statistics data. 
Associated J/?/? -|- 7 production The energy spectrum of J/ip produced in as- 
sociation with a photon, 7^ J/ip + '-fX, is a distinctive probe of colour-octet 
processes | 146|J127| , 147 , |148|| . In the colour-singlet channel and at leading-order in as, 
J/tp + J can only be produced through resolved photon interactions. The corre- 
sponding energy distribution is thus peaked at low values of z. The intermediate and 
large- 2; region of the energy spectrum is expected to be dominated by the colour-octet 
process •yg cc[8,^Si] 7. Colour- singlet and colour-octet contributions to associated 
J/ip + 'y photoproduction through direct and resolved photon processes are shown in 
Figure |l^ as a function of the J/ip energy. Observation of a substantial fraction of 
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Figure 17: Colour-singlet (CS) and colour-octet (CO) contributions due to direct and 
resolved photons to the J ftp energy distribution da/dz in associated J/tp + j production 
at the photon-proton centre-of-mass energy ^s^p = 100 GeV. The shaded area bounded 
by the solid lines represents the sum of all contributions with colour-octet matrix elements 
according to Equation 0. The separate colour-octet channel is plotted for {0'^/'^[8,^ Sq]) = 
(C)-^/V'[8,3Po])Mc = 0.012 GeV^. Other parameters as in Figure [ll. 

J/tp + J events at z ^ 0.5 would provide clear evidence for the presence of colour- 
octet processes in quarkonium photoproduction. For a more comprehensive study 
of associated J/ip + 'j production including various other differential distributions 
see ID-481 . 
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Bottomonium production With tlie significant increase in statistics at the up- 
graded HERA colhder, it might be possible to study inelastic photoproduction of 
bottomonium states for the first time. The large value of the bottom quark mass 
makes the perturbative QCD predictions more reliable than for charm production, 
and the application of NRQCD should be on safer grounds for the bottomonium 
system where ~ 0.1. However, the production rates for T states are suppressed, 
compared with J/^/^, by a more than two orders of magnitude at HERA, a consequence 
of the smaller bottom electric charge and the phase space reduction by the large h 
mass. The next-to-leading order colour-singlet cross section in the inelastic region 



> 1 GeV and z < 0.8 is shown in Figure as a function of the photon-proton 
centre-of-mass energy. Several years of data taking at the upgraded HERA may be 
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Figure 18: Next-to-leading order colour-singlet prediction for inelastic T(IS') photopro- 
duction at HERA as a function of the photon-proton centre-of-mass energy ^s^p. Param- 



eters: MRST parton distribution functions [129|; factorisation and renormalisation scale 
/i = \/2m;, ||5^. The error band is obtained by varying the bottom quark mass and in 
the range 4.5 GeV < rric < 5 GeV and 0.1225 < as(M^) < 0.1175, respectively. 



required to obtain accurate experimental information on the inelastic T photopro- 
duction cross section. 
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5 Conclusion and Outlook 



The NRQCD factorisation approach is a systematic framework for inclusive quarko- 
nium production and provides a striking and elegant explanation for the large direct 
J /ip and ipC^S) cross sections measured in pp collisions at the Fermilab Tevatron. 
Still, more theoretical work and experimental information is needed to firmly es- 
tablish the applicability of NRQCD factorisation to charmonium production. The 
absence of transversely polarised J/ip and ^^(25*) at large transverse momentum, if 
confirmed with higher statistics data, could imply that the conventional NRQCD 
power counting rules have to be replaced by alternative schemes, or that higher-order 
QCD effects arc essential to describe the spin-depcndcncc of the charm cross section. 

The substantial increase in luminosity at the Tevatron Run 11 will allow to assess 
the role of colour-octet contributions in bottomonium production more reliably. The 
calculation of ncxt-to-lcading order QCD corrections and a systematic treatment of 
soft-gluon radiation arc, however, required to describe the cross section in the pt 
region probed by the experiments at the Tevatron. A comprehensive test of the 
NRQCD factorisation approach in the bottomonium sector will have to wait until 
the LHC starts operating. If the charmonium mass is indeed not large enough for 
a non-relativistic approach to be reliable, the onset of transverse T polarisation at 
Pt 3> Mr may become the most decisive test of the NRQCD factorisation approach. 

A global analysis of different production processes is needed to verify the univer- 
sality of long-distance matrix elements in NRQCD. The colour-octet matrix elements 
obtained from Tevatron data are in qualitative agreement with analyses of other pro- 
duction processes, like charmonium hadroproduction at fixed-target experiments and 
charmonium production in B decays. However, the sizeable theoretical uncertainties, 
mainly related to QCD effects in the short-distance cross sections, prevent a more 
quantitative test at present. 

The experimental data on charmonium photoproduction are adequately described 
by the colour-singlet channel, including next-to-lcading order corrections in a^. But 
there is a considerable amount of uncertainty in the normalisation of the theoretical 
prediction, and no constraints on the size of colour-octet processes can be deduced 
from the photoproduction analysis at present. The upgraded HERA coUider has a 
remarkable potential to test the theoretical framework more thoroughly. The analysis 
of resolved photon processes and charmonium polarisation, in particular, could be 
essential in assessing the importance of different quarkonium production mechanisms. 

The theoretical issues to be addressed in the future include a critical evaluation of 
the power counting rules and the calculation of higher-order QCD corrections. The 
importance of non-vanishing initial state transverse momentum has to be clarified 



52 



by analysing various high-energy QCD processes. Moreover, systematic effects in- 
herent in NQRCD, such as the inaccurate treatment of energy conservation in the 
hadronisation of heavy quark pairs and the mismatch between partonic and hadronic 
phase space need to be addressed quantitatively. Finally, reliable estimates of non- 
factorisable corrections and production mechanisms beyond NRQCD are required to 
improve the theoretical analyses. Theoretical progress and the wealth of new results 
from the future Tevatron and HERA experiments will lead to exciting and compre- 
hensive phenomenology and result in a more complete and rigorous understanding of 
quarkonium physics. 
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